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ABSTRACT 

We compare our latest single and binary stellar model results from the Cambridge 
STARS code to several sets of observations. We examine four stellar population ra- 
tios, the number of blue to red supergiants, the number of Wolf-Rayet stars to O 
supergiants, the number of red supergiants to Wolf-Rayet stars and the relative num- 
ber of Wolf-Rayet subtypes, WC to WN stars. These four ratios provide a quantitative 
measure of nuclear burning lifetimes and the importance of mass loss during various 
stages of the stars' lifetimes. In addition we compare our models to the relative rate 
of type Ib/c to type II supernovae to measure the amount of mass lost over the entire 
lives of all stars. We find reasonable agreement between the observationally inferred 
values and our predicted values by mixing single and binary star populations. However 
there is evidence that extra mass loss is required to improve the agreement further, to 
reduce the number of red supergiants and increase the number of Wolf-Rayet stars. 
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> ] 1 INTRODUCTION. 

• I— I , 

^^ . There are many problems outstanding in the field of stel- 

$_( ' lar evolution. The most massive stars, those that end their 
5t I lives as a supernova (SN), the main questions are 'At what 
rate do stars lose mass?', 'Which stars give rise to which 
SNe?', 'What physical processes occur in binary systems?' 
and 'How important are rotation and magnetic fields?' One 
way to gain insight into these questions is to observe pop- 
ulations of stars and SNe. Massive stars are classified in a 
range of possible spectroscopic stellar types. The main cate- 
gories are blue supergiants (BSGs), red supergiants (RSGs) 
and Wolf-Rayet ( WR) stars. If we count the numbers of such 
stars in a single system (e.g. a stellar cluster or galaxy) we 
can test the accuracy of our stellar models by comparing to 
model predictions. 

When such comparisons are performed the agree- 
ment can be poor. For example the ratio of the BSG to 
RSG pop ulations has been an outstanding problem for 
some time (iLanger fc Maededll995l : iMaeder fc Mevnetll200ll : 



subtypes WN and WC stars a re predicted b y modern stella r 
evolution (|Mevnet fc Maedeiii2005 : Eldridge fc Vink|[200^ ') . 
Massive single stars end their lives as core-collapse SNe 
but the relation between the star and SN type is not straight 
forward. A list of core-collap se SN types includes IIP, IIL, 
lib, Iln, Ilpec, lb and Iqj (jTuratto. Benetti fc Pastorelld 
[2007]). The diflFerence between type II and type I SNe is the 
presence or absence of hydrogen in the SN spectrum. It is 
easy to identify which stellar models contain hydroge n or not 
but the subtype classifications are fairly arbitrary IIPopov 



iMaasev fc OlsenlbOOal ). The ratio of the number of RSGs to 
the number of WR stars is also not correctly predicted by 
stellar models l|Massey||2003l ). However the ratios of WR star 
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1993!; iDray fc T out 2003: Hcgcr ct al. 20031: lEldridge fc Tout 
.2004ij : lYoung. Smith fc Johnsoa .20051 '). The most direct 
method to link supernovae to stellar end point models is 
to searc h for the pr ogenitor star s of SNe in pre-explosion 
images l|Van Dvk et al. 2003; Sm artt el al.|[2004l . 

If we only split the SN types into II and Ib/c another 
comparison can be made between models and observations. 
The relative rates of these two SN types and how they 
vary with meta llicity have been e stimated observationally 
(|Prantzos fc Boissier 2003). iMevnet fc Maeder (2005) have 
shown that standard non-rotating single star models do not 



^ Type la SNe are thought to be thermonuclear events from 
explosive carbon burning in a degenerate carbon-oxygen white 
dwarf and are not considered in this work. 
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agree with these observed rates and that mass loss must be 
enhanced, in their case by rotation, to find an agreement 
between predicted SN ratios and those inferred from obser- 
vations. 

The studies mentioned so far only compare single 
star models to the observations. This ignores the observed 
fact that many stars are in binaries. Binary evolution has 



been studied widely , (e.g. Podsiadlowski. Joss fc Hsu 1993: 



200J: 



^, (e 



1 — r 



Wellstein fc Langerl 119 99: Belczvnski, Kaloeera fc 



Izzard. Ramirez-Ruiz fc Touti 



Buhkl 



2007 



Vanbeveren. Van Bever fc BelkusI 120071 : IVzquez et al.l 
Dray fc Tout! 120071 ). The methods employed to estimate 
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the effect of binaries is based either on rapid population 
synthesis or a small number of detailed models. Population 
synthesis employs formulae or table s based on detailed 
mode ls to follow the evolution, (e.g. iHurlev. Pols fc ToutI 
|200(J and the evolution of a star can be calculated in a 
fraction of a second rather than many minutes taken by 
a detailed code. The speed of calculations means that 
the full range of parameters that govern the outcome of 
binary evolution can be studied and their relative effects 
evaluated. The disadvantage is that in the complex phases 
of evolution, such as when the envelope is close to being 
stripped, the results of population synthesis can be quite 
spurious. 

We have constructed a set of binary stellar models cov- 
ering a wide range of bi nary parameter space us ing a detailed 
stellar evolution code ()Eldridge fc Toutll2004al '). With these 
models we have investigated the effect binary evolution has 
on the lifetimes of the various phases of stellar evolution 
and therefore the effect it has on the relative populations 
of massive stars and the relative numbers of different types 
of SNe. While the binary parameter-space resolution is low 
compared to rapid population synthesis studies it is still one 
of the highest resolution studies undertaken with a detailed 
evolution code. 

We note that we do not consider here the effect of rota- 
tion or magnetic fields. The effect of rotation on the struc- 
ture and evolution of stars is a com plicated p rocess and 
has been studied by Heger. Langer fc Wooslevi l|2000h and 
iHeger fc Lang cr (2000) and the extensive study beginning 
in lMevnet fc Maedcr. (.19971') and c ontinuing up to most re- 
centlv in iMevnet fc Maedeij (|2007l ). The effect of magnetic 
fields on stellar structure and evolution has also been in- 



yestig a ted bvlMevnet fc Maeder ( 2003) , Mevnet fc Maeder 
l|2004l '). lHeger. Wooslev fc Spruid (|2005l ). Our code does not 
include rotation therefore our study reveals the importance 
of duplicity alone. We discuss how rotation may affect our 
results in Section 6. 

In this paper we start by providing details of our stel- 
lar model construction. We then specify our definitions of 
the different phases of stellar evolution and discuss how bi- 
nary evolution affects the relative lifetimes of stars in various 
phases of stellar evolution. Using these lifetimes we predict 
how the BSG to RSG, the WR to O supergiant, the RSG 
to WR and WC to WN ratios vary with metallicity. We 
then predict how the relative rates of type Ib/c and type II 
SNe vary with metallicity. With each of these predictions we 
compare the values to the ratios inferred from observations. 
Finally we discuss the implications of our results. 



2 CONSTRUCTION AND TESTING OF THE 
STELLAR MODELS. 

We use a detailed stellar evolution code to construct sets 
of single and binary star models rather than a rapid pop- 
ulation synthesis code, because we require accurate treat- 
ment of complex stages of evolution, such as when the en- 
velope is close to being stripped, to obtain more accurate 
pre-SN models. We model the binary parameters with a 
reasonable resolution owing to the great amount of com- 
puter power that is now available. We have compared our bi- 
nary mo dels to the predictions of a rapid population synthe- 
sis code (llzzard fc Tout 20031 : llzzard. Ramirez-Ruiz fc Touti 
I2OO4I : llzzard fc Toutll2005l ) and found reasonable agreement 
for most of our predictions. However difficulty in assigning 
WR types to population synthesis helium stars meant the 
agreement was poor for the WC/WN ratio. 



2.1 Detailed Single Star Models 

Our detailed stellar models were calculated with the Cam- 
bridge stellar evolution code, STARS, originally developed 
by Egglcton ( 1 971f) and updat e d most recently bv lPols et al.l 
(j 19951 and lEldridge fc ToutI l|2004al '). Further details can 
be found at the code's web pageqj. Our models are avail- 
able from the same location for downloa d without restric- 
tion. They are similar to those described bv lEldridge fc Touti 
(2004bi) but we use 46 rather than 21 zero-age main-sequence 
models. Every integer mass from 5 to 4OM0 is modelled 
along with 45, 50, 55, 60, 70, 80, 100 and 12OM0. The 
initially uniform composition \s X = 0.75 — 2.5^ and 
Y — 0.25 -1- 1.5Z, where X is the mass fraction of hydro- 
gen, Y that of helium and Z the initial metallicity taking 
the values 0.001, 0.004, 0.008, 0.02 and 0.04. The compo- 
sition mixture is scaled solar and Z© — 0.02. Our models 
end after core carbon-burning or at the onset of core neon 
ignition. During these burning stages, the envelope is only 
affected to a small degree because tnuciear ^ ithermai. Fur- 
thermore the late stages of evolution are rapid and have a 
negli gible effect on the final observab le state of the progen- 
itor (| Wooslev. Heger fc Weave3l2002l ). 

Our opacity table s include the la t est low temperature 
molecular opacities of iFerguson et al.l l|2005h . This update 
leads to tiny changes in the radii and surface temperatures 
of red supergiant models because the altered opacity values 
only slightly modify the surface boundary conditions. 

All the models employ our standard mass-loss pre- 
scription becaus e it agrees best with observations of 
SN progenitors (lEldridge fc ToutI l2004bl) Th e prescrip- 
tion is based upon that of iDrav fc ToutI (|2003) but mod- 
ified in several ways. For pre-WR mass loss we use the 
rates of Ide Jager. Nieuwenhuiizen fc van der Huchtl (|l988l ) 
except for OB stars for which we us e the theoretical 
rates of IVink. de Koter fc Lamerg l|200ll) . When the star 
becomes a WR star (Xsu Tface < 0.4, lo g( reff/K ) > 
4.0) we use the rates of iNugis fc Lamerd (120001 ). We 
scale all rates with the standard factor (Z/Z b)°'^ 
(|Kudritzki. Pauldrach fc Pulslll987l : iHeger et al]|2003l ). ex- 
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cept for the rates of lVink. de Koter fc Lamera l|200H ) which 
include their own metallicity scaling. 

The most important change is that we scale the 
WR mass-loss rates with the by metallicity, eve n though 
this is not normally included. IVanbeverenI (|200l|) first in- 
cluded the scal i ng in a population synthesis model while 
ICrowther et al.l l|2002r ) first suggested that the scaling should 
be included from observations of WR stars in the LMC 
and in the Galaxy. Recent theoreti cal predictions also sug - 
gest the scaUng should be included (|Vink fc de KotenboOSh . 
Inclusion of this scaling results in a greater agreement 
with the change of WR population ratios with metallicity 
IjEldridge fc Vinkll2006l ). There is some uncertainty in the 
exact magnitude of the exponent we should use. To remain 
consistent we use the same exponent as for the non-WR 
stars. 

We synthesise a single star population by assuming a 
constant star formation rate and th e initia l mass function 
(IMF) of iKroupa. Tout and Gilmord (|l993l l. 



2.2 Detailed Binary Models 

Binary stars experience different evolution from that of sin- 
gle stars if their components interact. In a wide enough orbit 
a star is not affected by a companion. Duplicity allows the 
possibility of enhanced mass loss, mass transfer and other 
binary specific interactions (e.g. irradiation, colliding winds, 
surface contact, gravitational distortion) and hence a greater 
scope for interesting evolutionary scenarios. 

We have modified our stellar evolution code to model 
binary evolution. The details of our binary interaction algo- 
ri thm are relatively simple compared to the scheme outlined 
in lHurlev. Tout fc Pols (2002 ). We use their scheme as a ba- 
sis but we change some details which cannot be directly 
applied to our detailed stellar evolution calculation. We also 
make a number of assumptions in producing our code to 
keep it relatively simple. Our aim is to investigate the effect 
of enhanced mass-loss due to binary interactions on stellar 
lifetimes and populations, therefore we concentrate on this 
rather than every possible physical process which would add 
more uncertainty to our model. We also make assumptions 
in calculating our synthetic population to avoid calculating 
a large number of models. For example we do not model the 
accretion on to the secondary in the detailed code. We take 
the final mass of the secondary at the end of the primary 
code as the initial mass of the secondary when we create our 
detailed secondary model. This avoids calculating ten times 
more secondary models than primary models. 



2.2.1 Evolution of a binary 

In our binary models we always treat the primary as the ini- 
tially more massive star when and we only evolve one star at 
a time with our detailed code. When we evolve the primary 
in detail it has a shorter evolutionary timescale than the 
secondary which remains on the main sequence until after 
the primary completes its evolution and so we can determine 
the state of t he secondary using t h e sing le stellar evolution 
equations of iHurlev. Pols fc Toud (|2000t ) . When we evolve 
the secondary in detail we assume its companion is the com- 
pact remnant of the primary (a white dwarf, neutron star or 



black hole) and treat this as a point mass. We describe how 
we combine these models to synthesise a population below. 
To model a binary we begin by specifying a primary 
mass. Ml, a secondary mass, M2, and an initial separa- 
tion, a. We assume the orbit is circular and conserve an- 
gular momentum unless mass is lost from the system. We 
follow the rotation of the stars, considering them to ro- 
tate as solid bodies but we do not include processes such 
as rot at ionally- enhanced mass loss or mixing. We give the 
st ars initial rotation ra t es link ed to their initial masses as 
in iHurlev. Pols fc ToutI (|2000l ). When the stars lose mass 
in their stellar wind angular momentum is lost both from 
their spin and orbit. We ignore tides so stellar rotation and 
orbital rotation evolve independently until Roche-lobe over- 
fiow (RLOF) when we force the stars into synchronous ro- 
tation with the orbit, transferring angular momentum from 
the spin to the orbit or vice versa. 



2.2.2 Roche-lobe overflow 

In a binary system the equipotential surfaces are not spher- 
ical but are distorted by the gravity of the companion star 
and by the rotation of the system. If the stars are small com- 
pared to the size of the orbital separation these surfaces are 
nearly spherical but as the star grows relative to the sepa- 
ration the shape is distorted and becomes more ellipsoidal. 
A star's surface eventually reaches the LI Lagrange point 
where the gravity of both stars cancels exactly. If a star 
expands beyond this point then material begins to flow to- 
wards the other star. To include this in our models we deflne 
a Roche lobe radius, i?Li, such that the sphere of this radius 
has the same volume as the material within the Roche Lobe 
deflned by the equipotential surface passing through the LI 
point. Therefore when the star has a radius greater than the 
Roche lobe radius it transfers material on to the other star. 
We use the Roche lobe radius given bv lEggletonl l|l983h . 



0.49^1 



2/3 



(1) 



R-Ll ^ 

where gi = M1/M2 and a is the orbital separation. It is 
accurate to within 2% for the range < gi < co. When 
Ri > Rli we calculate the rate at which mass is lost from 
the primary according to. 



MiR = F{Mi)[ln{Ri/Ru)fMQyr-\ 

where 

F{Mi) = 3 X 10"^[min(Ml,5.0)]^ 



(2) 



(3) 



chosen by experiment to ensure mass transfer is steady 
l|Hurlev. Tout fc Polj[200^ ). 

Mass lost from the primary is transferred to the sec- 
ondary but not all is necessarily accreted. Accretion causes 
the star to expand owing to increased total mass and there- 
fore an increased energy production rate if M2 ^ M2/TKH 
where tkh is the thermal, or Kelvin-Helmholtz, timescale. 
We assume that the star's maximum accretion rate is de- 
termined by its current mass and its thermal timescale. 
We define a maximum accretion rate for a star such that 
M2 ^ M2/TKH. If the accretion rate is greater than this 
then any additional mass and its orbital angular momentum 
are lost from the system. In general stars with lower masses 
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have longer thermal timescales than more massive stars. Ef- 
ficient transfer is only possible if the two stars are of nearly 
equal mass so the thermal timescales are similar. This is an 
approximate treatment but provides a similar result to the 
more complex model of IPetrovic. Langer fc van der Huchd 
l|2005h who included rotation and found it led to inef- 
ficient mass-transfer. For compact companions we derive 
the maximum accretion rate from the Eddington limit 
l|Cameron fc Mocbil967i '). 



where 5a is the change in orbital separation and 5M\ is the 
mass lost in the time period considered, one model timestep. 
In this formalism 5M\ is negative so the orbit shrinks. By 
determining the mass-loss rate from the RLOF equation we 
find that the CEE ends naturally with the stellar radius 
shrinking once most of the envelope has been removed. 



2.2.3 Common- envelope evolution 

If RLOF occurs but does not arrest the expansion of the 
mass losing star, growth continues until the secondary is 
engulfed in the envelope of the primary. This is common en- 
velope evolution (CEE). During CEE it is thought that the 
envelope is ejected by some dynamical process and that the 
principal sourc e of energy for the envelope ejection is the 
orbital energy (|Paczvnskil [19761 : iLivio fc SokeiJll984 ). The 
orbital separation decreases and there is a chance that the 
two stars may coalesce before the envelope is ejected. Alter- 
natively the two stars are left in a close orbit, commonly one 
helium star for a massive progenitor and one main-sequence 
star. If the secondary then evolves to a helium giant a sec- 
ond CEE phase can occur possibly leading to a merger or a 
very close binary. 

CEE is modelled in population synthesis by first cal- 
culating the energy required for the ejection of the enve- 
lope, the binding energy, -Ei,bind = GMiA4'i,cnv/(A-Ri, where 
Mi^env is the mass of the envelope, Ri is the radius of the 
star and A is a constant to reflect the structure of the en- 
velope. This is then compared to the initial orbital energy, 
^-orb.i = —GMc,iM2 / {2ai) , where Mc,\ is the mass of the 
core and Oi is the initial separation at the onset of CEE. If 
there is not enough orbital energy to eject the envelope and 
leave a stable binary then the two stars merge to form one 
star. Otherwise t he envelope is removed lea ving a stable bi- 
nary system fe.g. lHurlev. Tout fc Polsll2002l ) . The new orbit 
is calculated by the equation i?i,bind = otcY,{Eorh,i — Eorb.f), 
where E^rh.f is the final orbital energy from which the final 
separation is calculated. There is some uncertainty in the 
constant ctcE as there are other sources of energy available 
such as the energy from nuclear reactions and re-ionisation 
energy of the hydrogen in the envelope so it can be greater 
than one. 

In our models when CEE occurs we derive the mass- 
loss rate from the above equation for RLOF, however we 
limit the mass-loss rates to a maximum of 10"'^ M0yr~^ be- 
cause more rapid mass loss causes the evolution code to 
break down. We assume that the secondary accretes no 
mass because the CEE occurs on the thermal timescale 
of the primary which is normally shorter than the ther- 
mal timescale of the less massive secondary. The accre- 
tion rate for a compact remnant is re stricted to the Ed- 
dington limit IjCameron fc Mockl 119671 ). To calculate the 
change in the orbital separation we calculate the binding 
energy of the material lost in the CEE wind ((5-Ebinding = 
— (Ml -\- M2)SM\/ R\) and remove it from the orbital energy 
(^„,bit = -GMiM2/(2o)). We find. 



Sa — 



a^ Ml + Kh 5Mi 



Ri Ma 



Ml 



(4) 



2.2.4 Binary mergers 

In the most compact binaries the stars can merge during 
CEE to form a single star whose mass is the sum of its par- 
ents. We find that some of our binary models enter CEE 
and the binary separation tends to zero. When this occurs 
we use a different binary model. The primary is evolved in 
the detailed code and once its radius is equal to the binary 
separation all the mass of the secondary is accreted on to 
the surface of the primary star and the evolution contin- 
ues as a more massive single sta r. This is an extremely sim- 
f le mo del. Other models such as lPodsiadlowski. Joss fc Hsul 
1993) show that a more detailed treatment may be required 
but the efficiency of the merger is always uncertain. Our 
models provide an estimate of the effects of stellar merg- 
ers and are an upper limit as we assume all the secon- 
daries mass is accreted. Furthermore the outcome depends 
on the evolutionary state of the primary, it is possible that 
fresh hydrogen could be mixed into the helium burning zone 
and therefore cause the envelope to be explosively removed 
(|Podsiadlowski. Morris fc Ivanovall2006r ). 



2.2.5 The binary model subsets 

We calculate three subsets of binary models which we then 
combine to simulate a population of stars. The first subset 
consists of our primary models where we evolve the primary 
star in the detailed code and evolve the secondary according 
to the single-star evolution equation of iHurlev. Pols fc Toud 
(|2000t ). We choose our grid to be Mi =5, 6, 7, 8, 9, 10, 11, 12, 
13, 15, 20, 25, 30, 40, 60, 80, 100 and 120; <? = ^ = 0.1, 0.3, 
0.5, 0.7 and 0.9; log(a/RQ) =1.0, 1.25, 1.5, 1.75,^2, 2.25, 2.5, 
2.75, 3, 3.25, 3.25, 3.5, 3.75 and 4;. We do not model lower 
masses because these stars are never luminous enough to 
affect our population ratios and do not lead to core-collapse 
SNe. 

The second subset consists of our secondary models 
where we evolve the secondary in the detailed code after the 
end of the primary evolution and replace the primary by 
its remnant. We chose three masses for the compact com- 
panion to model a white dwarf, neutron star or black hole, 
Mi,post-SN = 0.6, 1.4 or 2M0 respectively. We choose the 
grid to use the same separation and mass grid as for our pri- 
mary model subset. The final subset consists of our merger 
models and we choose this grid to have the same Mi, a 
and q distribution as our primary model. We only calculate 
models with initial orbital separation smaller than the max- 
imum stellar radius for each initial primary mass. These are 
only used if the orbital separation reaches zero in a primary 
model. 
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2.2.6 Synthesising the population 

The binary population is synthesised by first assuming 
the primary masses are distributed with the IMF of 
iKroupa. Tout and Gilmord (119931 ) . Then we assume flat dis- 
tribution over 1 < Iog(o/R0) ^ 4 and 0.1 ^ q ^ 0.9. We 
calculate the primary stars' lifetimes and SN types using the 
primary models. If a primary model during evolution enters 
a GEE and the binary separation shrinks to zero we replace 
the primary model with a merger model. 

The end point of the primary model is taken to be at 
the onset of neon burning, or the onset of thermal pulses in 
an AGB star. First we determine the remnant the primary 
will form. A white dwarf is selected unless the primary un- 
derwent a SN if there is an oxygen-neon core and the carbon- 
oxygen core mass is greater than 1.3 M0 when a neutron star 
is selected. A black hole remnant replaces this if the helium 
core mass is greater than 8 M0 . In the case of a white dwarf 
the system remains bound. When a neutron star is formed 
we consider the effect o f its natal kick on the binary orbi t. We 



Table 1. Our definitions of stellar types. 



use the prescription of iBrandt &: Podsiadlowski |l995l ) and 
the kick distribution of [Hansen fc PhinnevI ( 19971 ) to deter- 



mine whether the binary remains bound or not. In the black 
hole case we assume the kick velocity is z ero and estimate 
the rem nant mass and ejected mass as in lEldridge fc Touti 
l|2004bl ). 

If the system is unbound we calculate the contribution 
of the secondary to the stellar lifetimes and SN type from 
our single star models. The contribution from binaries that 
remain bound are taken from our secondary models. The ini- 
tial mass taken for the secondary star is the initial secondary 
mass unless the star accreted material from the primary in 
which case we use the new mass after mass-transfer. 

Gombining the model subsets creates our synthetic pop- 
ulation of binary stars based on detailed models from which 
we can estimate the relative numbers of different stellar 
types and the relative rates of different types of SNe. We 
calculate these populations at flve different metallicities, 
Z = 0.001, 0.004, 0.008, 0.02 and 0.04. 



2. 2. 7 Uncertainties 

Uncertainties arise from assumptions we have made in the 
physical model for binary systems. The greatest uncertainty 
is that we do not model the primary and secondary in detail 
simultaneously. This means we are not able to model sys- 
tems where q is close to one. In these mass transfer can be 
quite efficient with the seco ndary accreting a larg e fraction 
of the mass of the primary (|Cantiello et al.ll2007^ . Our ap- 
proximation is reasonable because for primaries with initial 
masses < 6OM0, the secondary stars with q < 0.9 remain 
on the main-sequence during the primary lifetime. The pri- 
maries > 6OM0 with companions that evolve off the main 
sequence before the primary experiences a SN also rarely 
interact as both stars lose the hydrogen envelopes before 
becoming RSGs. The error introduced by this approxima- 
tion is dwarfed by the uncertainty in our assumed initial 
parameter distribution for the binaries. If we only consider 
binaries closer than \og{a/ Rq) < 3 then we find that our 
RSG/WR ratio decreases by at least a factor of two. 
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Figure 1. Hertzsprung-Russel diagram indicating regions that 
we include as BSGs, RSGs and WR stars with single star tracks 
to indicate the mass ranges they sample. The numbers indicate 
the initial stellar mass in M0 for each track. The initial metalhc- 
ity was Z=0.02. The BSG and WR regions overlap in the lighter 
shaded region adjacent to the BSG region. The distinction be- 
tween WR and BSGs is made by surface hydrogen abundance. 



2.3 Defining Stellar Types 

Stellar types may be defined by observable characteristics 
such as colour and luminosity. However for a complete and 
accurate classification spectra are required. For example 
some objects may initially appear to be BSGs but their 
spectra indicate they are quite different and are in fact WR 
stars. WR spectra contain broad emission lines and weak or 
no hydrogen lines indicative of fast and dense stellar winds. 
In this work we use the theoretical stellar definitions that 
are summarised in Table [l] 

First, we only consider massive BSGs, RSGs and WR 
stars that have log(L/L0) ^ 4.9, equivalent to a luminosity 
cut-off of A/boi 1^ —7.5. If we compare our models to sin- 
gle clusters with well defined ages the limit would not be 
critical. However, because we are considering a less specific 
population with an assumed constant star-formation rate, 
lowering the limit slightly to \og{L/LQ) ^ 4.7 does not sig- 
nificantly change our predicted BSG to RSG and RSG to 
WR ratios but does decrease the WC to WN ratio because 
a larger number of low-mass helium stars are included as 
WR stars. 

BSGs and RSGs both have hydrogen envelopes but 
differ in their surface temperatures. BSGs are hot with 
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log(reff/K) ^ 3.9. This includes O, B and A stars. RSGs are 
cooler with log(Teff/K) < 3.66. This includes K and M stars. 
We show how these regions fit the tracks of our single star 
evolution models in Figure [T] RSGs are stars more massive 
than 12Mq while BSGs are stars more massive than 2QMq. 
Our luminosity limit re duces contamination from AGE stars 
(JMassev fc Olsenll2003l ') in the RSG sample. The limit is also 
greater than the luminosity of stars that experience a blue 
loop during helium burning. This feature is very sensitive 
to the details of conve ction included in the st ellar models 
and so best excluded (ILanger fc Maederlll995l ). Our mod- 
els with blue loops lie outside of our RSG definition so our 
predicted ratios are less sensitive to the complications intro- 
duced by the details of convection and mixing employed in 
our models. 

We also consider the O supergiants (OSGs) which are a 
subset of the BSGsfor the WR/0 ratio. OSGs are the hottest 
hydrogen-rich stars with log(reff/K) ^ 4.48 and a surface 
hydrogen mass fraction, X > 0.4. There is no lu minosity 
limit for these stars as in lMevnet fc Maederl l|2005l ). 

The other stellar types we consider are WR stars. We 
define a WR star to be any star that has log(reff /K) ^ 4.0 
and a surface hydrogen mass fraction, X < 0.4. We also 
require that \og{L/LQ) ^ 4.9. Further to this we subdivide 
WR stars into WNL, WNE and WG stars. In each case the 
element nitrogen or carbon is dominant in their emission 
spectrum. The sequence is due to the exposure of nuclear 
burning products on the surface of the WR stars. We assume 
a WR star is initially a WNL star. When X < 0.001 it 
becomes a WNE star. The star becomes a WG star when 
X < 0.001 and {xc + xo)/y > 0.03 where xc, xq and y are 
the surface number fractions of carbon, oxygen and helium. 



fecting the total stellar mass and thus the vigour of the 
nuclear reactions in the core. Secondly, stronger mass loss 
strips mass more quickly from the stellar surface, which leads 
to the exposure of hydrogen and then helium burning prod- 
ucts and results in the appearance of different WR subtypes. 

The WR lifetimes are longer at higher metallicities as 
the WR phase is entered at an earlier stage and there is 
a longer period of core helium burning during the WR 
phase because of larger pre-WR mass loss. The kink seen 
in the Z=0.02 plot (Fig. [2]) is due the fact that stars with 
masses below it undergo a red supergiant phase when mass 
is stripped and the surface hydrogen abundance drops below 
the limit set for WNL stars (Xsurfaco ^ 0.4). However, the 
stars are still cool and some time passes before they are hot- 
ter than the temperature limit (10* K). Stars more massive 
than the mass at the kink spend most of their evolution at 
temperatures hotter than this limit and become WR stars 
as soon as the surface abundance requirement is reached. 

Our primary model WR lifetimes in Figures [2] and [3] 
have relatively minor differences as far as the total WR life- 
times are concerned. The main difference is that the min- 
imum mass for a WR star is decrease to around 15Af0 at 
both metallicities. These stars will produce lower mass WR 
stars than are created by single stars. 

The lifetimes of the WR subtypes at Z = 0.02 indicate 
that the primary models have longer WNE lifetimes than 
single stars and shorter WG lifetimes. The WR stars below 
the minimum mass limit for single WR stars do not have a 
WG phase and exist only as WN stars. At Z = 0.004 we see 
the same trend but many of the lower mass WR stars only 
have WNL phases. 



3 BLUE AND RED SUPERGIANT AND 
WOLF-RAYET LIFETIMES. 

In this section we compare the stellar lifetimes of our sin- 
gle and primary models at two different metallicities. We 
calculate a mean lifetime from our primary models only to 
represent the effect on the lifetimes of RLOF and GEE. The 
mean lifetime for primary models with the same initial mass 
is calculated assuming a flat distribution in q and log a. 

The primary model BSG mean lifetimes are slightly in- 
creased relative to the single star BSG lifetime. The total av- 
erage lifetime for BSGs increases by 4 to 8 per cent and there 
is no trend with metallicity. This is because close merger sys- 
tems extend the burning lifetimes of the primary star. Also 
stars initially too low in mass to be counted as BSGs accrete 
enough mass to rise above our luminosity limit. 

The RSG lifetimes of our models are shown by the dot- 
ted lines in Figures [2] and |3l The single star lifetimes at 
Z = 0.004 are longer than those at ^ = 0.02. This is due 
to the reduction in mass-loss rates at the lower metallicity 
so more time is required for the hydrogen envelope to be re- 
moved. Our primary models have RSG lifetimes two to three 
times shorter. RLOF and GEE greatly reduce the number 
of RSGs by removing their hydrogen envelopes to form WR 
stars. 

The WR lifetimes of our s ingle stars are identical to 
those of lEldridge fc VinkI l|2006l ). Mass loss affects the life- 
times. Firstly, it determines the total WR lifespan by af- 



4 RATIOS OF STELLAR NUMBER COUNTS 

Gounting and comparing the number of stars in different 
stellar populations is a useful method for investigating stel- 
lar evolution. A comparison between observed populations 
and those predicted by a stellar evolution code is only worth- 
while if in the count stars are not missed or misclassifled. 
Many older observations relied on groups of stars deter- 
mined by photometry. The use of spectrometry has im- 
prove d the accuracy and co mpleteness of observational sur- 
veys (JMassev fc Olsenll2003l ). 

The metallicity mass fraction of the observed stellar 
populations was calculated from the log[0/II] + 12 values 
and comparing them to the values from our sets of models. 
Because this process is ambiguous and the position of solar 
metallicity for the mass-loss scaling is indistinct, we have 
assumed the metallicities are uncertain by 25 per cent. 



4.1 Blue to Red Supergiant Ratio 

The ratio of the number of BSGs to the number of RSGs 
has been a problem in stellar astrophysics for some time 
(|Langer fc Maedeijl 19951 ) . Observations tend to indicate that 
the BSG/RSG ratio should decrease with metallicity while 
most stel lar evolution codes predict a constant or increas- 
ing ratio. lEggenberger. Mevnet fc Maederl (|2002l ') suggested 
that the situation could be improved by using spectroscopy 
to confirm the identity of supergiants. If we compare our 
models to their observed trend we find good agreement 



The effect of massive binaries on stellar populations and supernova progenitors 7 



Single star models, Z = 0.02 



Single stars models, Z = 0.004 




40 60 80 

initio! Moss / Mq 

Primary models, Z = 0.02 



120 





20 40 60 80 100 

initio! Ivloss / IvIq 



Primary models, Z = 0.004 



120 




40 60 80 

initio! Ivloss / Ivlg 



20 40 60 80 100 

initio! Ivloss / Ivlg 



Figure 2. RSG, WNL, WNE and WC lifetimes at Z=0.02. Upper 
panel, our single star models. Lower panel, our primary models. 
The lifetimes of the binary stars are mean lifetimes assuming a 
flat q and log a distribution. Each region is labelled but progresses 
from bottom to top by WNL, WNE and WC phase. 



Table 2. Observed numbers and popula tion ratio of BSGs and 
RSGs. Values are taken from lMassev fc Olsen (2003 ). 



Figure 3. As for Figure[2]but Z=0.004. 



System 


log[0/H] + 12 


A^BSG 


-^^RSG 


A^bsgZ-'Vrsg 


SMC 
LMC 


8.13 
8.37 


1484 
3164 


90 
234 


16 ±6 
14 ±4 



at solar metallicity but no a Rreement at SMC metallicity. 
Ho wever Chiosi et al. (11993) indica t e tha t the SMC clus- 
ter lEggenberger. Mevnet fc Maederl (|2002l ) observed has a 
spread of ages similar to the age of the cluster itself, there- 
fore the ratio they observe is not due to a single population 
of s tars with the same age. 

iMassev fc QlsenI (|2003l ) provide a detailed observational 
study of RSGs in the Magellanic clouds with rigid definitions 
for BSGs and RSGs (see Table [T]). They found from observ- 
ing the entire SMC and LMC that the ratios of BSGs to 
RSGs were 16 and 14 respectively (see Table [2]). 

Comparing the observations to predicted values from 
our models in Figure |4] we find that single stars predict too 
low a value because either there are too many RSGs, too 



s 10 




0.004 0.008 

Metallicity 



0.04 



Figure 4. Ratio of the numbers of BSGs to RSGs versus m etal- 
licity. Observations are taken from lMassev fc Olsenl l l2003r . The 
solid line is from our single star models while the dashed line is 
from our binary models. 
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few BSGs or a combination of these two. Extra BSGs can 
be made via stellar mergers. Such stars are seen in globular 
clusters when the mass of the merged star is greater than 
the mass of the cluster turn-off, these are blue stragglers 
l|Sills et al.ll2000f ). The second binary effect is that stars be- 
low the minimum initial single star mass for a WR star will 
lose their hydrogen envelopes to become WR stars because 
of RLOF or GEE. These are stars that would have been 
RSGs and therefore the total number is reduced. In com- 
bination these processes increase the predicted BSG/RSG 
ratio as our binary models show. 

The binary model ratio is dependent on how we choose 
the distribution of binaries in q and log (a/R©). In our bi- 
nary population we include very wide binaries, log (a/R©) Js 
3, which evolve as single stars so our binary population 
is a mix of interacting binaries and single stars. Approx- 
imately one third of our primary models evolve as single 
stars. Therefore, to reproduce observations, two thirds of all 
stars must be in interacting binaries. If our binary models 
are not completely correct something extra, such as rotation, 
may decrease the r equir ed fraction of interacting binaries. 
iMaeder fc Mevnetl l|2001f ) find that rotation has a strong ef- 
fect on the BSG to RSG ratio but the trend they found was 
that the BSG to RSG ratio decreased rather than increased. 
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Figure 5. Ratio of the numbers of Wolf-Rayet stars to O- 
supergiants versus met allicity. Observations are taken from 



iMaeder fc Mevnetl 1 11994) . The solid line is our single star models 
while the dashed line is our binary mo dels. The dashed-dot t ed lin e 
for the Geneva models is taken from iMevnet fc Maedeii l|2005|) . 
The y-axis erro r bars are an ass umed error of 50 percent of the 
values given bv IMaeder fc Mevneti (^994 ). 



4.2 Wolf-Rayet to O supergiant Ratio 

The WR to OSG ratio has been studied for sometime 
IMaeder fc Mevnetl [l993 ') . If binaries leave the RSG pop- 
ulation to become WR stars therefore this ratio should in- 
crease if the BSG to RSG ratio decreases. The observed ratio 
is less certain than the BSG to RSG ratio because there is 
;reater uncert ainty in the completeness of the observations 
Massevll2003l ). However the trend that the ratio decreases 
with metallicity is well established (JCrowthcr 2007). 

Figure [S] shows that our single star models underpredict 
the ratio while our binary models are in better agreement, 
as are the Geneva rotating models. This is because in our 
models more stars are stripped of their envelopes and be- 
come WR stars. There is also a contribution to the OSG 
population from stellar mergers and secondary accretion as 
for the BSG/RSG ratio above. The trend our binary models 
predict with metallicity is a little too shallow but within the 
uncertainty of the observed ratios. This again could indicate 
something extra may need to be included in our models. 



4.3 Red Supergiant to Wolf-Rayet Ratio 

The RSG to WR ratio compares the relative populations 
of the two stellar types that have completed core hydrogen 
burning and so measures the influence of mass loss. 

In Table[3]we list the observed populations. We plot the 
ratio in Figure [6] it decreases dramatically with increasing 
metallicity. Such a change would require a stronger scaling 
of mass loss with metallicity than we currently employ. 

Our single star models agree with the observations at 
SMC metallicity while our binary models agree with the ob- 
servations around LMG metallicity. One way to match the 
observed trend between these metallicities is to have a metal- 
licity dependent binary fraction. However something extra is 
still required around solar metallicity to get an exact agree- 
ment. At the higher metallicities the ratios are based on a 



Table 3. Observed numbers and population ratios of RSG and 
WR stars. The values are taken from Massev (2003). 



System 


log[0/H] + 12 


A^RSG 


A^WR 


A''r,sg/'^wr 


SMC 


8.13 


90 


12 


7.5 ±4.7 


NGC6822 


8.25 


8 


4 


2.0 ±1.8 


M33(o) 


8.35 


26 


18 


1.44 ±0.95 


LMC 


8.37 


234 


108 


2.17 ±0.87 


M33(m) 


8.6 


7 


22 


0.32 ±0.25 


M33(i) 


8.75 


3 


17 


0.18 ±0.16 


M31 


9 


1 


19 


0.053 ±0.058 



small number of observed stars (see Table [Sj and therefore 
do not sample the full binary parameter space. Furthermore, 
small numbers mean our assumptions of constant star for- 
mation and IMF become invalid. For example, if all the WR 
stars are in very close binaries the RSG/WR ratio would be 
even smaller than we estimate. A large number of stars must 
be observed to calculate the population ratio to ensure the 
full range of possible binary systems is probed. 



Table 4. Observed numbers and population ratios of WC and 
WN stars. The values are taken from Mcynot fc Maeder (2005). 



System 


log[0/H] + 12 


A^wc 


AfwN 


iVwc/A^WN 


SMC 


8.13 


1 


11 


0.083 ±0.10 


NGC6822 


8.25 





4 


0.000 ±0.20 


M33(o) 


8.35 


3 


16 


0.190 ±0.17 


LMC 


8.37 


17 


91 


0.190 ±0.11 


M33(m) 


8.6 


12 


35 


0.370 ±0.25 


MW{< 3kpc) 


8.7 


30 


35 


0.880 ±0.52 


M33(i) 


8.75 


11 


19 


0.580 ±0.42 


M31 


9 


13 


14 


0.900 ±0.66 
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Figure 6. Ratio of tlie number of RSG to WR s t ars ve rsus metal- 
licity. Tlie observations are taken from iMassevI 1120031) . The solid 
line is from our single star models while the dashed line is from 
our binary models. 
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Figure 7. Ratio of the number of WC to WN stars versus metal- 
licity. The observations and the dashed-dot t ed lin e for the Geneva 
models are taken from iMevnet fc Maeden l l2005h . The solid line 
is from our single star models while the dashed line is from our 
binary models. 



4.4 WC to WN Ratio 

Finally, we consider a ratio which provides a measure of the 
relative lifetimes of the two main WR star types: WN and 
WC (we include WO with WC). The lifetimes of the two 
types are determined by the mass-loss rates of WR stars. 

We plot the observations and predicted ratios in Fig- 
ure [7] it shows our predicted ratios and the most re cent 
Geneva group rotating models IjMevnet fc Maedeiil2005l ) for 
comparison. The trend of the observations is for the ratio 
to decrease with m e tallici ty. In agreement with the results 
of lEldridge fc VinkI l|2006l ) we find that the single-star mod- 
els that scale the mass-loss rate of WR stars with initial 
metallicity agree with the observed trend. 

The binary models in this case give a lower WC/WN 
ratio than the single star models. This is because, as can be 
seen in Figures[2]and|3] the lifetimes of WN stars increase by 
a greater factor than the lifetimes of WC stars in binaries. 
If we combine a population of single and binary stars the 
resulting ratio is too low and requires an increase in the 



WC population relative to the WN population to regain an 
improved agreement. The result i s simil ar to that found by 
IVanbeveren. Van Bever fc Belkui l|2007h . 



5 SUPERNOVA PROGENITORS 

There are many studies which investigate the con- 
nection between SNe and massive stars. Some stud- 
ies consider single-star evolution and predict the ini- 
tial parameter spa ce and th e rela ti ve rate of dif- 
ferent SN types jHeger et all l2003l : lEldridge fc Toud 
l2004bl : lllirschi. Mevnet fc Maedej|2004) . while other stud- 
ies are concerned with the evolution of binary-stars 



JPodsiadlowski. Joss fc HsJl993 : de Ponder fc VanbevererJ 
120031 : llzzard. Ramirez- Ruiz fc Tou t 2004). In this section we 
use our single- and binary-star models to predict the rela- 
tive SN rates and determine how they vary with metallicity. 
Firstly, we link our models to each SN type then, secondly, 
we predict how the relative SN rates vary with metallicity. 

Core-collapse SNe are classified according to their 
lightcurve shapes and spectra. Matching stellar models to 
ob served SNe i s difficu lt an d example scheme s can be found 
in iHeger et all (|2003h and lEldridge fc Toud (|2004bh . Here 
we check the amount of hydrogen in the progenitor model: 
if there is more than O.OOIMq of hydrogen left in the stellar 
envelope the SN is of type II, otherwise type Ib/c. 

There are many subtypes of SNe. The main distinguish- 
ing criterion is the mass of the hydrogen or helium enve- 
lope at the time of explosion but sometimes the circum- 
stellar environment must als o be considered l|Heger et al.l 
l2003l : lEldridge fc Toudl2004bl '). In this paper we group SN 
II sub-types (e.g. P, L) together as type II and SN lb and 
Ic together as type Ib/c. In our single-star models the ini- 
tial mass rang e of Ib/c progenitors is r estricted to the most 
massive stars l|Eldridge fc Toudl2004bl ). In our binary-star 
models the full range of masses can lead to type Ib/c SNe. 

SN rates in different ga laxy types hav e been 
measured for some time |Cappellaro et al.l Il997l : 



ICappellaro. Evans fc Turattdll999l ) and, more recently, SNe 
observations have been used to determine how the relative 
rate of type Ib/c to type II SNe varies with metallicity 
l|Prantzos fc Boissied 120031 ). The errors especially in the 
absolute rates of the searches are considerable owing to the 
small sample size and the uncertainties in the completeness. 
The relative rates are less uncertain as the selection effects 
are of similar magnitude and cancel. 

We plot our predicted SN rate ratios against the ob- 
served ratios in Fig. |8] The observations indicate a general 
trend of a decreasing rate of type Ib/c SNe relative to type 
II as metallicity decreases. This is as expected owing to the 
decreasing strength of stellar winds with reduced metallic- 
ity meaning that fewer stars lose all the hydrogen before 
core-collapse. 

We find that our theoretical predictions agree with 
the trend indicated by the observations. However the sin- 
gle star models predict a lower relative rate for type Ib/c 
SNe than the observations and the binary star models pre- 
dict a value that agrees with the observations. We also plot 
the SN rate ratio predicted by the Geneva rotating mod- 
els (|Mevnet fc Maededl2004l ) which agree with observations. 
The conclusion we draw is that we must consider rotating 
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Figure 8. The observed and predicted ratios of the type Ib/c 
SN rate to the typ e II SN rate. Observations are taken from 



iPrantzos fc Boi ssien (12 0031). The Ge neva model predictions are 
taken from lMevnet fc Maedeii 1120051 ). the upper Hne is for their 
rotating models while the lower line is from their non-rotating 
models. 



models and/or binary star models to exactly simulate the 
observations. 



6 DISCUSSION &: CONCLUSIONS 

We have compared our stellar models of massive stars to sev- 
eral observations and find that by including binary stars we 
achieve better agreement. This is because duplicity increases 
the chance of a star losing its hydrogen envelope which leads 
to fewer RSGs, more WR stars and more type Ib/c SNe than 
single star models predict as required by observations. 

Problems remain and the agreement is not perfect. The 
BSG/RSG ratio and RSG/WR ratio indicate that we stiU 
predict too many RSGs or not enough WR stars. Also the 
binary models predict to many WN stars and too few WC 
stars. Physics we have not included in our models might 
resolve these problems. The main culprit is rotation. It 
has two effects on our predictions. First it introduces mix- 
ing that extends the time stars spend on the main se- 
quence and therefore increases the number of BSGs. Sec- 
ond, if rotation is rapid it enhances mass loss by reduc- 
ing the depth of the potential from which mass must es- 
cape and therefore increases the population of WR stars 
at the expense of the RSG population. Rotationally en- 
hanced mass loss affects the WC/WN ratio by shortening 
the WN lifetimes. However the effect of rotation can be 
mu ch more complex then t hese simple effects. For exam- 
ple (MaedeE^iMe^ne^ llooj) find that the BSG to RSG ra- 
tio is decreased rather than increased by rotation owing to 
a change of the internal helium gradient at the hydrogen- 
burning shell. In addition tides in a binary may boost the 
importance of rotation during importa nt phases of evolution 
and complicate the situatio n further j[Petrovic et al.ll2005l : 
IPetrovic. Langer fc van der Hucht 2005|). 

Another limitation of our models is that we do not 
include mass losing eruptions that are common in the 
luminous-blue variable (LBV) stars. These most massive 
stars experience giant outbursts, losing large amounts of 



mass in one short event. T his has dramatic implic ations for 
the evolution of the object (jSmith fc OwockilboOa ). There is 
growing evidence that SN progenitors can experience these 
outbursts prior to th e SN explosion (Kotak fc Vink 2006|; 
IPastorello et al.ll2007h . Such outbursts would have a similar 
effect on stellar populations as rotation. They reduce the 
time stars spend as RSGs and could reduce the time taken 
to remove the helium envelopes for WN stars. Therefore a 
better understanding of these outbursts may also improve 
agreement between stellar models and observations. 

In summary, binary stars must be considered when com- 
paring stellar evolution models to observations. However, 
binary evolution alone cannot explain all the observations. 
Fine tuning of stellar models is still required and the effect of 
enhanced mass loss owing to rotation or LBV-like eruptions 
must be further considered. 
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